Temperature and Field Dependence of the Anisotropy of MgB 2 
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The anisotropy 7 of the superconducting state of high quality single-crystals of MgB2 was de- 
termined, using torque magnetometry with two different methods. The anisotropy of the upper 
critical field was found to be temperature dependent, decreasing from 7 ~ 6 at 15 K to 2.8 at 35 K. 
Reversible torque data near T c reveal a field dependent anisotropy, increasing nearly linearly from 
7 ~ 2 in zero field to 3.7 in lOkOe. The unusual temperature dependence is a true bulk property 
and can be explained by non-local effects of anisotropic pairing and/or the k— dependence of the 
effective mass tensor. 

PACS numbers: 74.60.Ec, 74.20.De, 74.25.Ha, 74.70.Ad 



The discovery of superconductivity at T c « 39 K in 
MgB2 |Q has caused a lot of interest into it's physi- 
cal properties (for a review see Ref. Q). Measurements 
of the isotope effect |jl and e.g. the n B nuclear spin- 
lattice relaxation rate §J indicated a BCS type s-wave 
phonon-mediated superconductivity. Calculations of the 
band structure and the phonon spectrum predict a dou- 
ble energy gap [|, ||, a larger gap attributed to two- 
dimensional p x — y orbitals and a smaller gap attributed to 
three-dimensional p z bonding and anti-bonding orbitals. 
A substantial body of evidence by, among others, spe- 
cific heat measurements 0, point-contact spectroscopy 
scanning tunneling spectroscopy ^| and penetration 
depth measurements [jl0| has emerged to support this 
scenario. An alternative scenario with a single, but 
anisotropic, gap was also proposed (TTJ , and supported 
by recent Raman measurements J^] . 

A double gap structure or an anisotropic energy gap 
should influence the anisotropy in the superconducting 
state, according to the standard anisotropic Ginzburg- 
Landau (GL) theory 7 = (m*/TO* fc ) 1/2 = A c /A ah = 
Cab/Cc = H^/H^i where ||a6(||c) indicates the field H 
perpendicular (parallel) to the c— axis of the sample and 
m* , A, £ and H C 2 are the GL effective mass, the penetra- 
tion depth, the coherence length, and the upper critical 
or bulk nucleation field, respectively. Most values re- 
ported for the anisotropy of polycrystalline or thin film 
MgB2 span the range of values of 7 = 1.1 — 3 B, but 
there are also reports with 7 rs 6 — 9 |l3|, [lij. Up to 
now, there are only four reports on transport measure- 
ments of the upper critical field anisotropy performed on 
single crystals, giving values of 2.6 p5[ , 2.7 jl6| and 3 
Jl7| , ft8| . Magnetic measurements of the angular depen- 
dence of H C 2(0), yielding 7 = 1.6, were reported only on 
aligned crystallites Q . 

Here, we report magnetic torque measurements on sin- 
gle crystals of MgB2, performed in a wide range of tem- 
peratures from 15 K to 36 K in magnetic fields of up to 
90kOe. We provide evidence that the bulk anisotropy 



7 is not universally constant, but is temperature depen- 
dent down to at least 0.4 T c and shows a pronounced field 
dependence near T c . Microscopic origins of the unusual 
T— dependence of 7 in MgB2 are discussed. 

We have grown single crystals of MgB2 with a high 
pressure cubic anvil technique, similar to the one de- 
scribed in Ref. jl6j . The details of crystal growth will be 
published elsewhere. In brief, a mixture of Mg and B was 
put into a BN container and a pressure of 30 — 35 kbar 
was applied. Growth runs consisted of heating during 
lh up to the maximum temperature of 1700 — 1800°C, 
keeping the temperature for 1 — 3 h and then cooling to 
room temperature during 1 — 2 h. Flat crystals were up 
to 0.8 x 0.6 x 0.05 mm 3 in size, with sharp transitions to 
the superconducting state at about 38 — 39 K. 

Measurements were performed on miniaturized 
piezoresistive cantilevers specifically designed for torque 
magnetometry [po[ . The torque f=rnxB^mxH, 
where rn is the magnetic moment of the sample, was 
recorded as a function of the angle between the applied 
field H and the c— axis of the crystal for various fixed 
temperatures and fields. For measurements close to T c , 
in fields up to 14 kOe, a non-commercial magnetometer 
with very high sensitivity was used. For part of these 
measurements, a vortex-shaking process was employed 
to speed up the relaxation of the vortex lattice plj . The 
observation of a well-resolved lock- in effect in t(9) (see 
upper inset of Fig. |5|) indicates there are no variations 
of crystallographic alignment throughout the samples. 
A crystal with a volume of about 4 x 10~ 4 mm 3 (sample 
A) was measured in this system. Another crystal with a 
volume of about 8 x 10 -3 mm 3 (sample B) was measured 
in a wider range of temperatures down to 15 K in a 
Quantum Design PPMS with torque option and a 
maximum field of 90 kOe. 

An example of a torque vs. angle curve is given in the 
inset of Fig. [I]. For small angles 9 the torque is essentially 
zero. Only when H is nearly parallel to the ab— plane 
there is an appreciable torque signal. The curve can be 
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FIG. 1: Torque r vs. angle data, showing the definition of 
8 C 2- The inset shows the full torque curve measured. 



interpreted in a straight-forward way: for H parallel to 
the c— axis the sample is in the normal state, while for 
H parallel to the ab— plane it is in the superconducting 
state. The crossover angle 9 C 2 between the normal and 
the superconducting state is the angle for which the fixed 
applied field is the upper critical field. The inset of Fig. 
[l] also shows hysteretic behaviour due to irreversibility. 
The irreversibility field H m (T,9) determined from the 
torque measurements is very high, close to H c i- Prelimi- 
nary SQUID measurements on similar crystals indicate a 
much lower H m J22]j ; an extended discussion of the irre- 
versible properties of MgP>2 will be published elsewhere. 

The crossing between straight lines through the back- 
ground and the superconducting torque signal was used 
to define C 2. This definition is not unambiguous. Taking 
the analysis of the data more strict we have to apply the 
appropriate scaling rules. The magnetization M of a 3D 
system in the GL theory of fluctuations in the vicinity of 
the transition temperature T c (H) in high magnetic fields 
is given by a universal function F of the distance from 
T C (H) @: 



M _ T 2 / 3 / A(T - T c ) 
~H ~ WI* \ (TH) 2 / 3 



(1) 



where A is a material constant. Combining the above 
dependence with the angular dependence of the torque 
p3 we find that the rescaled torque signal 



P = -re 1/3 {9) / (sin^ cos 9H 5/3 (l - l/-/ 2 )T 



2/3 



(2) 



with e(9) — (cos 2 9 + sin 2 0/7 2 ) 1 / 2 , is a universal func- 
tion of the distance from T c with a fixed value F(0) at 
T = T C (H). Taking into account the F(0) value for the 
theoretical dependence of the universal function for a 3D 
system p5[ we can estimate that for a volume of the sam- 
ple of 8 x 10 _3 mm 3 P reaches at T — T C {H) a value of 
about 2 x Kr 10 dyncmOe~ 5/3 KT 2/3 . The inset in Fig. 
^ presents the angular dependence of the rescaled torque 
P in different magnetic fields at 22K. The crossing of the 
P(9) dependence for each field with the line of the con- 
stant value of 2 x 10~ 10 dyncmOe _5/,3 K _2 ^ 3 determines 
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FIG. 2: Upper critical field H c2 vs. at 22 K. The full line is a 
free fit of Eq. ([|) to the data. Alternative fits with fixed values 
of the anisotropy 7 are also shown. The inset shows selected 
rescaled torque P (see Eq. (^)) vs. angle 0. From left to right, 
curves shown were measured in H = 24, 30, 40, 55, 70, 80, and 
85 kOe. The two criteria used for the determination of C 2 are 
indicated by dashed lines. The data in the main panel were 
obtained employing the lower criterion. 



the H c2 (9) dependence as it is shown in the main panel of 
Fig. ||. It is important to stress that the results obtained 
depend not very sensitively on the criterion chosen and 
it will be shown later (see Fig. |J) that even with a three 
times higher criterion we get very similar temperature de- 
pendences of H C 2 and 7. Additional t(H) measurements 
at fixed angle give H C 2(9) values corresponding well to 
those from t(9) measurements. 

Within the applicability of the anisotropic GL theory 
the angle dependence of the upper critical field is pre- 



dicted to be 1 26 



H c2 (9) =H^ (cos 2 9 + sin' 



oh 2 ) 



-1/2 



(3) 



A fit of Eq. (|) to the data at 22 K yields 7 = 5.1(1) and 
ffj 2 c = 17.2(1) kOe. This fit (7 = 5.1) describes the data 
well, while alternative fits with 7 fixed to 4 and 6 are 
clearly incompatible with the data, as shown in Fig. ||. 
Figure shows the angular dependence of H c i scaled 

lie . 

by Hjg to directly compare the anisotropy at different 
temperatures. The 15 K data are well described by the 
line corresponding to 7 = 6, while the 34 K data lie below 
the line for 7 = 3.5. The data indicate an anisotropy sys- 
tematically decreasing with increasing temperature. To 
show this is not an artifact related to fitting, we present 
the angular dependence of the rescaled torque P for fixed 
H/H^2 m the inset. The curves clearly shift to higher 
angles with increasing temperature. Furthermore, we di- 
rectly checked t(9) raw data in fields above and below 
and b to give absolute limitations of 7. We find 
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FIG. 3: Angular dependence of H C 2(8)/H'^, employing the 
lower criterion, at various temperatures. Curves correspond 
to H C 2(6) according to Eq. (^|) for 7 = 3,3.5,4.5, and 6. In- 
set: Some rescaled torque P (see Eq. (Q)) vs. angle 8 curves, 
corresponding to H/H^ — 2.67. 



that at 22 K, the anisotropy must be higher than 4.4, 
while at 34 K, it must be lower than 3.5. 

All data are summarized in Fig. |4|. The data ob- 
tained from fits to Eq. (|3|) do not vary much with the 
criterion used for the determination of 9 C 2, and agree well 



with calculated from thermal conductivity data |2S|] 
measured on a single crystal grown with the same tech- 
nique. The T— dependence of H}% is in agreement with 
calculations by Helfand et al. [£9|. The corresponding 
H'j^O) ~ 31kOe is relatively small compared to litera- 
ture values, which may indicate that the crystals inves- 
tigated are relatively free of defects. The values 

lie 

were obtained from the two fit parameters H^ 2 and 7. 
There is a slight positive curvature of ij]'^ h (T), which 
can be attributed to the T— dependence of 7. The 7(T) 
dependence may also be the origin of the positive cur- 
vature of H C 2 observed in other measurements of bulk, 
thin film and single crystal measurements B]. Due to 
the lack of low temperature data and the variation of 
7, only a rough estimation ifj.'^O) w 230 kOe can be 
given. The anisotropy data show that 7 systematically 
decreases with increasing temperature. A change of the 
criterion used for the determination of 9 C 2 leads to small 
shifts of the magnitude of 7, but the temperature de- 
pendence is always the same. The highest upper critical 
field anisotropy 7 ~ 6 was obtained at 15 K, the lowest 
anisotropy 7 ~ 2.8 at 35 K. From Fig. || we estimate 
7 (0) = 7 - 8 and 7 (T C ) = 2.3 - 2.7. 

Small systematic deviations from Eq. (|3|), observed 
near T c , indicate that the field influence on 7 may be 



important as well. To clarify this point, we have mea- 
sured the reversible torque r as a function of angle 9 for 
various fields and temperatures near T c . The data were 
analyzed with an equation derived by Kogan et al. , 
based on the anisotropic London model, which contains 
the GL anisotropy 7 as a parameter. 

Without shaking, the irreversibility of the torque was 
relatively large and a clear lock-in effect was observed. 
The inset of Fig. || shows an example of shaked torque 
data, together with the fitted curve. An evaluation of the 
data measured for various T and H up to lOkOe with 
the Kogan formula |^7j reveal that 7 is field-dependent 
with larger 7 in larger fields, while temperature varia- 
tions do not affect 7 appreciably (see Fig. |^). Indications 
of a field dependence of the coherence length £ (a precon- 
dition of a field dependence of it's anisotropy) have been 
observed by specific heat (3^] and muon spin rotation |31 
measurements on NbSe2- We should stress that Kogan's 
formula assumes equivalence of the anisotropies of the 
coherence length £ and the penetration depth A. Since 
the anisotropy of A might be different from the one of £ 
in our case, the two different methods used in this work 
can lead to different effective values of 7. The field de- 
pendence of 7 may be related to the peculiar double gap 
structure of MgB2 with a large gap of two-dimensional 
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FIG. 4: a) Upper critical field H C 2 vs. temperature T. Open 
symbols correspond to H \\ ab, full symbols to H \\ c, from fits 
of H C 2 (0) data to Eq. Q . Up triangles are from measurements 
on sample A (with 8 C 2 determined as shown in Fig. |l|) and 
squares (circles) are from measurements on sample B, using 
the lower (higher) criterion of constant rescaled torque P (see 
inset of Fig. 0). b) Temperature dependence of the upper 



critical field anisotropy H^ 2 
H c2 {6) to Eq. (|). 
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FIG. 5: Field dependence of the anisotropy 7 as determined 
from the angular dependence of the reversible torque t(9). 
Upper inset: raw data, demonstrating the effect of shaking 
and the lock- in effect. Lower inset: fitting of the shaked data 
to a formula derived by Kogan et al. M, yielding 7 ~ 3.2. 



the BEM tensor and the anisotropic EPC, and is beyond 
the scope of this work. 

In conclusion, the upper critical field anisotropy 7 of 
MgB2, determined by torque magnetometry, decreases 
with increasing temperature. Measurements of the re- 
versible torque near T c reveal an almost linear field 
dependence of the anisotropy of the coherence length 
and/or the penetration depth as well. Our results im- 
ply a breakdown of standard anisotropic GL theory with 
a (temperature and field independent) effective mass 
anisotropy. The temperature dependence of 7 can be 
approximated tentatively by j(T) = 7* + 7(1 — T/T c ) n 
with n close to 1. Here, 7* ~ 2.3 — 2.7 is the band ef- 
fective mass anisotropy and 7 ~ 4.5 — 5.5 arises from the 
anisotropy of the attractive electron-electron interaction 
and/or the wave- vector dependence of the effective mass 
tensor. 

We thank B. Batlogg, P. Miranovic, A. Sologubenko, 
and I. L. Landau for useful discussions. This work was 
supported by the Swiss National Science Foundation, 
by the European Community (contract ICAl-CT-2000- 
70018) and by the Polish State Committee for Scientific 
Research (5 P03B 12421). 



nature and a small three-dimensional gap, which is very 
rapidly suppressed in a magnetic field g. 

A temperature dependent iJ c2 anisotropy was previ- 
ously observed, e.g. in NbSe 2 f| and LuNi 2 B 2 C [||. 
However, in MgB 2 , the effect is much more pronounced. 
It was shown that any theory capable of explaining a tem- 
perature dependence of 7, needs to take into account non- 
local effects HI, [35| , which can be pronounced in samples 
of high purity. In the vicinity of T c , non- locality is not 
important Therefore, m* c /m* ab = 7 2 (T C ) ~ 5-7 cor- 
responds to the standard GL effective mass anisotropy. 
We note that this is significantly higher than the calcu- 
lated |5|, ^6) anisotropy of the band effective mass (BEM) 
averaged over the Fermi surface (1.0 — 1.2). The micro- 
scopic theories |34|, show that the T— dependence of 7 
cannot be attributed to an anisotropy of the BEM tensor, 
unless it is also wave vector dependent. An anisotropic 
energy gap, caused e.g. by an anisotropy in the electron- 
phonon coupling (EPC), can also lead to variations of 7 
with temperature. 

In MgB 2 , first principles calculations suggest both a 
pronounced wave-vector dependence of the BEM tensor 
p6| and a strong anisotropy of the EPC (see e.g. || [37j). 
The latter leads naturally || to the observed double en- 
ergy gap and was suggested Q to be responsible for 
the unusually high T c of MgB 2 . To our knowledge, there 
is only one theory calculating a temperature dependent 
H c2 anisotropy of MgB 2 pT| , which predicts, however, 

/ to increase with increasing T. A quantitative 
explanation of the measured 7(T) apparently needs to 
take into account both the wave-vector dependence of 
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